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' Almost since the time it was formulated, the overwhelming consensus has been that random matrix theory 

, (RMT) is in excellent agreement with neutron resonance data. However, over the past few years, we have 

CO . obtained new neutron-width data at Oak Ridge and Los Alamos National Laboratories that are in stark dis- 

agreement with this theory. We also have reanalyzed neutron widths in the most famous data set, the nuclear 
data ensemble (NDE), and found that it is seriously flawed, and, when analyzed carefully, excludes RMT 
with high confidence. More recently, we carefully examined energy spacings for these same resonances 
in the NDE using the A3 statistic. We conclude that the data can be found to either confirm or refute the 
theory depending on which nuclides and whether known or suspected p-wave resonances are included in the 
analysis, in essence confirming results of our neutron-width analysis of the NDE. We also have examined 
radiation widths resulting from our Oak Ridge and Los Alamos measurements, and find that in some cases 
they do not agree with RMT. Although these disagreements presently are not understood, they could have 
broad impact on basic and applied nuclear physics, from nuclear astrophysics to nuclear criticality safety. 



^ . 

O ■ 1 Introduction 
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In 1956, Porter and Thomas (T\ proposed a theory to explain the surprising discovery that s-wave resonance 
reduced neutron widths (Fj^) spanned a very wide range. Starting from seemingly sound and fundamental 
assumptions, their theory predicted that r° values follow a distribution with one degree of freedom 
{Vn = 1), which subsequently became known as the Porter- Thomas distribution (PTD). 

Almost since that time, the overwhelming consensus has been that data and theory agree very well. In 
fact, faith in this theory is so strong that in the past ~30 years it has been extremely rare to find a paper in 
the literature in which new data have been used to test the theory. Instead, standard procedure has been to 
use the theory to correct new data for experimental deficiencies. In the intervening years, random matrix 
theory (RMT) Q was developed and has placed the PTD on more formal footing, broadened the scope 
and predictions, and provided links between nuclear physics and many other fields, including quantum 
chaos. As a consequence, the impact of neutron resonance data has become much broader, as such data are 
routinely cited as some of the best proof of the veracity of RMT. 

In this paper, we summarize recent tests of RMT predictions using neutron resonance data. Because all 
experiments from which such data are obtained have important limitations, we begin by discussing these 
problems and how they are circumvented in Section |2] 

Over the past few years, we have obtained new Fj'^ data |l3]|4l at the Oak Ridge Electron Linear Acceler- 
ator (ORELA) 15J and Los Alamos Neutron Science Center (LANSCE) ||6l facilities that are much better 
than previous data. We discuss results of using these data to test the PTD in Section[3] 
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The 72^-matrix analyses from which FJ; data are obtained also yield total radiation widths (F^) for these 
same resonances. These data also can be used to test the PTD and other predictions of theory, but the 
situation is complicated by limited knowledge of the level density as a function of excitation energy and 
the photon strength functions needed to model the 7 decay of the resonant capturing states. In Section]?] 
we discuss some general observations regarding what has been learned from data obtained at ORELA 
over the past several years, as well as a surprising very recent result |7:| discerned from our ^''^Sm data 
taken at LANSCE. 

In addition to our own data, we have reanalyzed |]8] neutron widths in the most famous data set ll9l [T0l . 
the so-called nuclear data ensemble (NDE). More recently, we have reanalyzed energy spacings for these 
same resonances, with careful attention to evaluation of uncertainties. We discuss results of these analyses 
in Section|5] 

In almost all cases we have examined so far, we find significant deviations from the PTD or other 
disagreements with RMT predictions. There also have been several other reported disagreements with the 
PTD ifTTHTSll . which we briefly describe in Section |6] To our knowledge, those results have never been 
explained or retracted. 

As far as we know, these disagreements with RMT presently are not understood, although our recent 
results have (re)ignited theoretical interest IIT6HT91 , some of which is closely related to earlier work iBOl 
on non-Hermitian Hamiltonians. We very briefly discuss recent theoretical work in Section |7] Because 
many models assume that both particle and primary gamma transitions follow the PTD, violation of this 
assumption could have broad impact on basic and appUed nuclear physics, from nuclear astrophysics to 
nuclear criticality safety. 

Reliable separation of small s- from large p-wave resonances remains the largest barrier to better tests 
of RMT using neutron resonance data. In Section |8] we very briefly describe new experiment techniques 
aimed at surmounting this barrier. 

We end with a short summary and some conclusions in Section |9] 

2 Experiment limitations 

An ideal test of the PTD or other predictions of RMT would involve a complete (no missing resonances) 
and pure (all resonances having the same parity) data set. Typically, s-wave neutron resonances are used 
because they are largest and hence easiest to observe. Unfortunately, all experiments have a lower limit for 
observing resonances as well as a threshold below which s- and p-wave resonances cannot be differentiated. 
These facts are illustrated in Fig. [1] in which neutron widths for i^^Pt resonances from our ORELA 
measurements are plotted as a function of resonance energy. 

To understand how the data in this figure illustrate these facts, it is useful to know the following details. 
Because ^^^Pi is a zero-spin target (I'^ = 0+), all s-wave resonance have spin and parity J'" = i^, and 
statistical spin factor gj = (2/+i)(2j+i) ~ ^' where j = ^ is the neutron spin. In contrast, there are two 
spin/parity values, ^ and | , for p-wave resonances, for which gj = 1 and 2, respectively. Shown in 
Fig. [T]are effective reduced neutron widths gF° = (?F„/ y/En/leV. For s-wave resonances, these are 
the usual reduced neutron widths which should, on average, be constant. However, for p-wave resonances, 
these effective reduced widths should be, on average, proportional to En- Also, according to standard 
level-density models, the average number of resonances should be proportional to 2 J + 1 and independent 
of parity. Therefore, it is expected that there should three times as many p- as s-wave resonances (i.e., 
Dq = 3Z)i). Finally, due to the overall larger penetrability for s waves and the fact that Pt is near the peak 
of the s- and valley of the p-wave neutron strength function, s-wave resonances should be much larger than 
p-wave ones at these energies. 

These expectations are in qualitative agreement with the data in Fig. [T] For example, it is evident that 
there is one group of resonances with larger and constant (on average) gF^ and another, more populous 
group, having smaller effective reduced neutron widths with sizes (on average) proportional to En- We 
were able to make firm s-wave assignments to almost all resonances in the former group by virtue of their 
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asymmetric shape (due to interference with the comparatively large s-wave potential background) in the 
total cross section. Such resonances are depicted by solid symbols in this figure. Also, we were able to 
make a few firm — ^ (and hence p-wave) assignments under the assumption that the total 7 width 
does not vary by more than a factor of three. These resonances are marked by X's in Fig. [T] Unfortunately, 
both techniques for assigning resonance parity cease to work when the neutron width becomes too small. 
Nevertheless, it is obvious from this figure that if an analysis threshold shown by the solid curve is used, 
there is only an extremely small chance that a p-wave resonance will be included. However, it is still likely 
that some s-wave resonances are below this threshold, so this fact must be accounted for in any subsequent 
analysis. These data represent some of the best ever acquired for the purpose of testing RMT; by virtue of 
the excellent sensitivity for observing small resonances, the good separation of s- from p-wave resonances, 
and the number of firm parity assignments. Even so, it is vital to realize that a pure and complete set of 
s-wave resonances still could not be obtained. 

To avoid systematic errors, these experiment facts must be incorporated into the analysis technique. 
These problems affect r°, F-,, and resonance-spacing data in different ways, so each will be discussed 
separately below. 

3 Testing the PTD with r° data 

The maximum-likelihood (ML) technique is most often used to test whether data are consistent with the 
PTD because it is statistically efficient. Because the PTD is a special case (;/„ — 1) of the family of 
distributions, it is assumed that reduced neutron widths are distributed accordingly and the ML method is 
used to estimate the most likely value of Given the shape of the distribution as a function of i^n, 
neglecting the effect of missed s-wave resonances below threshold will, in general, lead to a falsely large 
value of t'n from the ML analysis. Conversely, including even just a few p-wave resonances in an s-wave 
set will, in general, lead to a falsely small value of !^„. 

Typically, these difficulties have been surmounted by using an energy-independent threshold as an in- 
tegral part of the ML analysis, implicitly assuming that all s-wave resonances above threshold have been 
observed. We recently have shown Q that an energy-dependent threshold (on r° ) of the form T — aEn, 
where a is a constant factor and £"„ the neutron energy, offers three advantages compared to using an 




5000 10000 15000 



Fig. 1 Effective reduced neutron widths versus energy for Pt resonances from our ORELA data. Open circles, 
filled circles, and X's depict data for all, firm s-, and firm p-wave resonances, respectively. The threshold used in the 
ML analysis of Ref. @j is shown as a solid curve. 
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energy-independent threshold. First, p-wave contamination is eUminated equally effectively at all ener- 
gies. This is because the penetrability factor for p waves differs from s waves by (to good approximation) 
a factor of En- Second, experiment thresholds have approximately this same energy dependence; thus 
possible diffusiveness of the instrumental threshold can be surmounted equally effectively at all energies. 
Third, statistical precision of the analysis is maximized by allowing the largest p-wave-free set of s-wave 
resonances to be included. 

The analysis technique was described in Refs. |l4|[8l. Each resonance A has an energy E\ and reduced 
neutron width F"^. The probability density function (PDF) f{x\vn) for a distribution is given by: 

f{x\vn)dx ^ ^J^'L^^ (~i~) ' cxp(^-^)da;, (1) 

where G(^) is the gamma function for x — > F^,j/(F"), and (F°) is the average reduced neutron 
width. 

The joint PDF for statistical variables F^^^ and E\ is defined in a 2D region defined by inequalities 
Ex < i^max and F^^^ > T{E\), where E'max is an upper limit of energies Ex. The expression for this PDF 
reads 



^ An 



(2) 



The factor C, ensuring a unit norm of h^, is Vn- and (F") -dependent. The ML function was calculated 
from all rig pairs [i?^'''^, T^''^] within the specified region obtained from the experiment. Specifically, 

no 

L {^n, K)) = n ^° [K"^ I (r° )) . (3) 



One problem with many previous tests of the PTD is that they were made using data from nuclides near 
the peaks of the p- and valleys of the s-wave neutron strength functions. As a result, neutron widths for 
the two parities are nearly equal and hence it is very difficult, if not impossible, to obtain a pure s-wave 
set spanning a sufficiently wide range of neutron widths to be useful for testing the PTD. For this reason, 
it is much better to work with nuclides near the peaks of the s- and valleys of the p-wave neutron strength 
functions. We recently have made measurements on several such nuclides; ^^^Sm and I92,i94,i95,i96pj 

Our ^'*^Sm(n, 7) measurements |3] were made using the Detector for Advanced Neutron Capture Ex- 
periments (DANCE) ll2n at LANSCE. One problem working with odd-A nuclides is that there are two 
possible s-wave spin states ( J = 3 and 4 in the case of ^^^Sm), and performing the best test of the PTD 
requires knowing all the resonance spins. Using 7-ray multiplicity data obtained with DANCE, we devel- 
oped O] and improved ll22l an new method for separating the two s-wave spins, and hence determined 
spins for nearly all observed resonances for En < 1 keV. This allowed us to use the combined J = 3 and 4 
data to test the PTD. As shown in Fig. |2] the surprising result was that the data changed from agreeing very 
well with the PTD (i/„ = 0.91 ± 0.32) for resonances below 350 eV, to significant disagreement with the 
PTD (I'n = 3.19 ± 0.83) for resonances within the next 350 eV. This change occurs at the same energy as 
a previously observed ll23i non-statistical behavior in the a-particle strength function ratio in this nuclide. 

Q^j. I92,i94,i95,i96pj ^j^j^ g| ^vgj-g taken at ORELA. Neutron capture and total cross sections for 
isotopically-enriched samples were measured using a pair of CeDg detectors and a ®Li-loaded glass de- 
tector, respectively. The neutron total cross section for a "°*Pt sample also was measured. Resonance 
parameters were determined from a simultaneous 7?.-matrix analysis of the nine sets of capture and total 
cross sections. Parameters for a total of 1264 resonances were determined, 631 of which could be assigned 
as s wave, accounting for nearly all resonances above the thresholds used in the ML analyses. The ^^^Pt 
data were excluded from the main ML analysis because spins could not be determined for enough reso- 
nances. The i92,i94pj. j^j^ each exclude the PTD by nearly three standard deviations (i^„ = 0.57±0.16 and 
0.47 ± 0.19, respectively). Results for i^^Pt were consistent with i92,i94pj.^ albeit with reduced statistical 
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precision — 0.60 ± 0.28) due to the smaller number of resonances observed for this isotope. Combined 
in a very conservative manner, the i92,i94pj. ^^^^ exclude the PTD with at least 99.997% confidence. 

4 Testing the PTD with data 

The total radiation width for a neutron resonance is the sum of partial widths r^-^ for the different 
channels by which the capturing state can decay by 7 emission, T-^ = X]"=i ^ii- According to RMT, 
partial widths F;^ follow the PTD (i^i^ = 1). In contrast to the neutron case where there is only a single 
open neutron channel (elastic scattering) at the relevant energies, there are many open channels through 
which 7 decay can proceed. Due to the properties of distributions, in the simplest model total radiation 
widths are predicted to follow a distribution with degrees of freedom given by the effective number 
of independently-contribution channels, rics = v-^ ~ ^"=1 ^«7- "^cff ~ 100, Y-^ distributions are 
predicted to be very narrow. Moving beyond this simplest theory requires detailed simulation of the partial 
radiation widths for primary transitions, which in turn requires accurate knowledge of the relevant level 
densities and photon strength functions. On the other hand, because distributions are much narrower 
than the PTD, it might be easier to detect deviations from theory, especially given the limited number of 
resonances typically available. 

Comparisons (2511 we have made to date using our ORELA data for Pt isotopes indicate that measured 
F-y distributions are broader than those predicted using standard photon-strength-function and level-density 
models. In addition, as shown in in Fig. [3] distributions for several nuclides measured at ORELA appear 
to contain extra tails. Although we have observed that simulations within the framework of the nuclear 
statistical model can lead to significant differences from x^ distributions predicted by the simplest model, 
we have been unable to simulate tails as large as those shown in Fig. [3] Hence, these tails may be the sign 
of non-statistical or collective effects. 

Also, as we very recently reported in Ref. Q and illustrated in Fig. 31 both the average total radiation 
width (F^) and the width of the F^ distribution increase rather abruptly near En = 300 eV for resonances 
in ^*^Sm. The median (variance) test 1261 indicates the null hypothesis that medians (variances) of the 
F^ distributions in the two energy regions are the same can be rejected at the 99.8% (99.9%) confidence 
level. In addition, the Smirnov and Cramer-von Mises two-sample tests 1261 reveal the null hypothesis that 




Fig. 2 Cumulative r„ distributions from our LANSCE data for "^Sm resonances for two energy regions. Shown are 
the number of resonances having reduced neutron widths larger than some value, versus that value. Circles and X's 
depict data for _B„ < 350 and 350 < En < 700 eV, respectively. Solid and dashed curves depict the PTD and a 
distribution for = 3.19, respectively. See text and Ref. |[3] for details. 
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data in the two energy regions were sampled from the same distribution can be rejected with > 99% and 
> 99.9% confidence, respectively. In essence, all these statistical tests indicate that the change in the 
distribution evident in Fig. |4]is highly statistically significant. 

Theoretical interpretation of this change may be aided by estimation of distribution parameters for the 
two regions. To this end, we used the ML method. As noted above, in the simplest model data are 
expected to follow a distribution with many degrees of freedom, typically i^-, ^ 100. For such large 
values of ly^, a distribution is very close to Gaussian in shape. One advantage of using a Gaussian rather 
than distribution for the analysis is that data uncertainties AF^ can easily be included jZTl . 

Therefore, we used the ML technique described in Ref. lIZTl to estimate most likely values for the 
means (F^) and standard deviations (Tn of the F^ distributions in the two energy regions. Resulting ML 
estimates are cttv = 4.67 ± 0.81, (F^) = 52.0 ± 1.1, and cttv = 11-7 ± 1.5, (F-^) = 59.6 ± 2.0, for 
the lower- and upper-energy regions, respectively. Hence, these ML results, which take into account the 
measurement uncertainties, also indicate that F-y distributions in the two energy regions are significantly 
different. Translated to distributions, these ML results yield i/^ = 248 ± 87 and 52 ± 14 for the F.y 
distributions in the lower and upper energy regions, respectively. 

These changes in the ^"'^Sm F^ distribution occur at the same energy that the neutron-width distribution 
changes for resonances in this nuclide, as described in Section[3]above. Furthermore, as shown in Fig. |5] 
these changes in the F° and F-^ distributions appear to be mirrored by increases in the both the average 
^"'^Sm(n, 7) cross section and fluctuations around the average near this same energy. Explaining these 
energy-dependent effects represents a significant challenge to theory. 



5 Testing RMT with the NDE 

The NDE ||9l[l0l is a set of 157 proton and 1250 neutron resonance energies consisting of 30 sequences in 
27 different medium- weight and heavy nuclides. The ensemble was assembled to test predictions of RMT. 
Fluctuation properties of resonance energies in the NDE were found to be in remarkably close agreement 
with RMT predictions for the Gaussian orthogonal ensemble (GOE). Although there have been several 
other successful tests of RMT using nuclear resonances, the NDE is perhaps the most important because, 
as stated in Ref. Q, "As a result of these analyses, it became generally accepted that proton and neutron 
resonances in medium weight and heavy nuclei agree with GOE predictions." Hence the NDE routinely 




Fig. 3 Cumulative F-, distributions from our ORELA data for three different nuclides. Shown are the fraction of 
resonances with F^,/ (F^,) larger than a given value versus the value, where (F^,) is the average F^, determined from 
the data. Solid staircase plots depict measured data whereas dashed smooth curves show distributions. 
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Fig. 4 Cumulative distributions for neutron reso- 
nances in ^''^Sm. Shown are the fraction of resonances 
with larger than a given value versus the value. Stair- 
case plots depict the measured data whereas smooth 
curves show Gaussian distributions from ML analyses. 
See text and Ref. (T) for details. 
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Fig. 5 Open circles depict our DANCE ^''^Sm(n,7) 
cross sections averaged over 80-eV-wide bins. Error 
bars corresponding to one-standard-deviation statistical 
uncertainties are smaller than the symbols. The solid 
and dashed curves show results of statistical model cal- 
culations based on the average resonance parameters for 
the lower- and upper-energy regions, respectively. See 
Ref. (7) for details. 



is cited as providing striking confirmation of RMT. We have carefully analyzed both neutron-width and 
neutron resonance-energy data in the NDE and found significant disagreements between these data and 
RMT predictions. These analyses, results, and an overall assessment of the NDE are described in the next 
three subsections. 

5 . 1 Testing the PTD with the NDE 

RMT for the GOE predicts that T'^ data follow the PTD. Reduced neutron widths have been reported for a 
subset of 1245 resonances in the NDE of Ref. ifTOl , consisting of 14 to 178 measurements for 24 nuclides. 
Details of the analyses, nuclides included, primary references, number of resonances in each, etc. can be 
found in Ref. d. 

Analyzing a data set comprised of many different nuclides such as the NDE involves at least two addi- 
tional potential pitfalls to those mentioned in Section |2] above. First, apparent sensitivities of the various 
experiments from which the NDE was derived differ by several orders of magnitude. Therefore, if the 
entire NDE were analyzed as a single set (as was done in Ref. IfTOll ). to avoid systematic error the threshold 
must be at least as high as the highest apparent individual threshold. However, doing this will exclude so 
many resonances from the analysis that the statistical precision of the result will be substantially reduced, 
and hence at least partially negate the reason for assembling the NDE in the first place. Therefore, separate 
thresholds should be used for each of the NDE nuclides. Second, the average reduced neutron width (r° ) 
depends on the shape of the distribution, and this shape may be different for each NDE nuclide. Therefore, 
it is important to include separate (F" ) parameters for each NDE nuclide. For these reasons, separate ML 
analyses must be made for each NDE nuclide, and then the combined result subsequently compared to 
theory. 

For the initial analyses, thresholds just below the smallest observed resonance for each NDE nuclide 
were used. Because experiment thresholds might not be precisely sharp, it is expected that the resulting 
i^n values would be systematically a bit large. However, almost all f„ values were less than the PTD 
value of 1.0. The weighted average of results at minima thresholds for the 24 nuclides in the NDE is 
i^n — 0.801 ± 0.052, which is 3.8 standard deviations smaller than the predicted result of i^n = 1- Hence 
these data reject the PTD with a statistical significance of 99.98%. A value significantly less than the 
PTD could be a sign of interesting physics. However, a more likely explanation is that the NDE contains 
sizable p-wave contamination. 
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That the NDE is contaminated by p-wave resonances is evident in many ways. First, many resonances 
in the NDE have been identified as p-wave or of uncertain parity. For example, in Ref. Il28l . 58 p-wave 
resonances in ^'^■^Th were assigned on the basis of 7-cascade information, 13 of which are in the NDE. One 
of these 13 resonances also is known II29II30I to be p-wave by its observed parity- violating asymmetry. For 
over half (14/24) of the NDE nuclides we analyzed, these resonances account for 5% or more of the total, 
for 10 of the 24 they are at least 10%, and in the three worst cases about 35%. Second, as explained in 
Ref. m, that many of the NDE resonances are, in fact, p wave is reinforced by the behavior of the Vn values 
from the ML analyses as functions of threshold. In many cases, Vn systematically increases with threshold 
before gradually stabilizing. This is just the behavior expected for a population of s-wave resonances 
contaminated by p-wave resonances. 

As explained in Ref. HI, removing effects of these p-wave resonances from the NDE ML analyses is a 
simple matter of raising thresholds until they are above the largest previously identified p-wave resonance 
and/or Vn stabilizes as a function of threshold. The resulting weighted average for the NDE is still in 
conflict with the RMT prediction for the GOE, albeit in the opposite direction from the result using the 
lowest thresholds: Vn ~ 1.217 ± 0.092, corresponding to a confidence level of 98.17% for excluding the 
PTD. Hence, when the NDE is cleansed of p-wave resonances, the data still reject the PTD with high 
confidence. Furthermore, because our ^^^Sm and i92,i94pj j-ggQj^g show that f„ can be on either side of the 
PTD value of 1.0, it seems likely that any analysis from which an average Vn is calculated from results for 
several nuclides may underestimate deviations from the PTD. 

In the next subsection, we describe our analysis of the NDE using the A3 statistic. This analysis 
illustrates several more related problems in the NDE. 

5.2 Reexamining A3 for the NDE 

The initial claim of excellent agreement between the NDE and RMT was based mainly on the A3 resonance- 
spacing statistic. Because our analyses described above indicated that the neutron-width data in the NDE 
do not agree with the PTD, we also examined the A3 statistic for this subset of the NDE data. We will 
refer to this subset of 1246 neutron resonances as the original NDE, or set I. There are 1246 energies but 
only 1245 neutron widths because one resonance in ^^^^^ has a reported energy but no width. 

As mentioned above, several of these NDE resonances have been assigned as either definite or probable 
p-wave resonances. Specifically, definite p-wave assignments, in a primary reference and/or one of the 
two compilations 13111321 . include the following 35 resonance energies (in eV): 3333.0 (^^"'Cd), 1702.8 
(i^^Yb), 1634.6 (i^'^Yb), 73135, 100220 and 125130 (^''Zn), 196.13, 391.53, 400.86, 411.62, 420.92, 
476.30, 540.09, 573.46, 764.7, 820.9, 850.5, 1114.9, 1204.1, 1372.54, 1387.7, and 1848.6 C^^TYl), and 
263.91, 454.1, 555.9, 732.5, 778.8, 1028.6, 1131.1, 1298.1, 1316.5, 1532.3, 1565.1, 1795.5, and 2070.9 
(-238^-) 'pjjg original NDE after exclusion of above-specified 35 resonance will be referred to as a corrected 
NDE, or subset Icorr- 

5.2.1 A3 statistic 

Given a selected interval of neutron energies {Ea,Eb), the A3 statistic is conventionally defined as 



Here, n{E) is equal to the number of the resonances in the interval {Ea,Ei,), while L = {Ei, — Ea)/E[D], 
where E[D] is an expectation value of spacing D between energies of neighboring resonances. 

Practical implementation of the statistic according to Eq. ^ faces the problem that, except for the 
maximum value of L for a given set of data, there is no unique definition of A3. That is, except for the 
maximum value of L, there are at least two pairs of Ea and Eb for which A3 can be calculated. To over- 
come this difficulty, we used a modified statistic, referred hereafter to as A3. For z-th sequence consisting 




(4) 
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of Ki energies Ef'^ < e!^'^ < ■ • • < i^A- of consecutive resonances of a given nucleus we introduced 
a function A3*''(L|Ao) of integer argument L > 2 and integer parameter Aq satisfying a condition that 
< Ao < - 1. Specifically, for 1 < L < - Ao + 1 



Af{L\Xo) = — ^ -ymin / [n,{E)-a - dE, (5) 

while for L < 1 or L > Aj — Aq + 1 the function A3'''(i|Ao) equals zero. 

Following the approach of Refs. ||9][l0l, all the individual A3 values for a given L were averaged over 
the range of possible Aq values to obtain an averaged statistic, (A3(L)), which is applicable to any subset 
I' of the full set X of the available experimental NDE sequences for individual nuclei; 

$:''^'''a;«(L|Ao) 
{A',{L)) = '^^'^7; — , (6) 

Quantities represented by right-hand sides of Eqs. (2) and dSj, as well as the average quantity defined 
by Eq. (|6]l, depend implicitly on resonance energies. In the spirit of the GOE model, these energies are 
assumed to be random variables, so that functions according to Eqs. ^ and (|6| are indeed statistics. 
Hereafter, depending on the context, we will refer them either to as functions or statistics. 

From a fixed subset I' C I of the NDE we took the resonance energies and determined from them func- 
tions ni{E) for all i e I' . Then, employing Eqs. (|5]l and (|6]l, we determined what we call an experimental 
realization of the averaged statistic {/S!^^"^ [L)) . 



5.2.2 Calculation of confidence limits 

Although the NDE is a fairly large data set, statistical uncertainties in {/S!^^^{L)) can still be substantial. 
To evaluate these uncertainties and assess whether these data are in agreement with theory, we employed a 
Monte Carlo technique to generate a large number (typically 5 x 10'^ to 10^) of artificial functions (A3 (L)) 
using energies E'j''' < i?2*'' < • • < i?]^'^ generated by the GOE model for all i G I'. In this way, we 
constructed an empirical probability distribution function for (A3(L)) at each value of L. We compared 
these simulations to the (A3'''^^(L)} values calculated from the NDE data to assess whether the NDE is 
compatible with predictions of the GOE model. 

Our implementation of the Monte Carlo technique used N x N GOE matrices. Hence, the N{N + l)/2 
independent elements of each matrix were statistically independent normally distributed random variables 
with zero mean. Variances of the off-diagonal elements were equal to a common value v, while those of 
the remaining elements were equal to 2v. 

(i) (i) (i) 

Generating energies E{ < E2 < • ■ • < Ej^ involved the following steps: 

1. For a given i E I' we chose a value Ni > max(2 A^, 100), adjusting for convenience the variance v 
to a value = l/4A^i. The largest value of Ni used was less than 1000. 

2. We created a random GOE matrix A'*^ of size Ni x Ni and determined its eigenvalues in ascending 
order, i.e. ef^ < e^'^ < •■• < e^|. 



3. Using the expression 



"27r 



2 arcsin (e^*^) + sin (2 arcsin(e^''' ) 



(7) 
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Monte Carlo simulations under the assumption that Fig. 7 Function {A'^"'''^ (L)) calculated from the cor- 

GOE model holds. Regions between pairs of thin curves rected NDE. See the caption of Fig. |6]for definitions of 

situated symmetrically to the median curve: ±lcr, ±2a the various plots. 



and ±3(7 confidence regions characterizing fluctuation 
of (A3(L)). Dashed line: function {A'^iL)) for the 
case of uncorrelated neutron resonance energies. 

where A = 1,2,..., A;, we transformed the set of eigenvalues {e^'''} to a modified set {^x^}- This 
transformation ensures that values of {ej'^} for iV^ — 5- cxd are distributed uniformly. 

4. Finally, from a part of values of {e^*''}, centered around zero, we created neutron energies. For this 
purpose we adopted a prescription E^^^ = ^\!,+\ ^'^^ A = 1, 2, ... , A^, where Aq satisfies a condition 
that e^*^^ < —Ki/(2Ni) < and that the length A^ is equal to that for the i-th sequence of the 

NDE of interest. 

In this way, we simulated a large number (typically 2x10'^ to 10^) of (A3(L)) values, for which the 
number and the lengths of the sequences are identical to those belonging to either subset I' C X or: 
I' C Icorr, and from the distributions of these values we calculated confidence limits for (A3(L)). 

Through further extensive Monte-Carlo simulations we verified that the limited sizes of the matrices 
used resulted in less than 0.2% systematic uncertainty in estimating confidence limits. 

5.2.3 Results 

Results of our analysis for case I' = I are shown in Fig. |6] where the function {/S!^^'^[L)) is plotted 
together with confidence limits of (A3(L)) . This figure indicates that the original NDE is in remarkably 
good agreement with GOE predictions for the entire allowed region L < 178, essentially confirming the 
results claimed in Ref . ||9l . 

The longest sequences of neutron resonance energies in the original NDE belong to ^'^^Th and ^'^^U. 
Results of our analysis limited to these two nuclides are illustrated in Fig. [8] Again, agreement between 
the data and GOE predictions is excellent. 

However, when ^^^Th and ^ssjj removed from the original NDE, the agreement disappears, as 
shown in Fig. |9] In this case, validity of the GOE model can be rejected with a confidence level greater 
than 99.99%. This conclusion is not at variance with the results shown in Figs. |6] and |8] but reflects the 
large weights with which ^■^^Th and ^^*U, by virtue of their comparatively long sequences, contribute to 
{A'^''^'^{L)). For example, for L > 70 this weight is greater than 80%. 



Copyright line will be provided by the publisher 




45 90 135 180 



L 

Fig. 8 Function {A'^"^''' (L)) calculated from the origi- 
nal NDE for the subset formed by ^^^Th and ^^*U res- 
onances. See the caption of Fig.[6]for definitions of the 
various plots. 
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Fig. 9 Function (A3™^(L)) deduced from the original 
NDE after eliminating the data belonging to ^'^^Th and 
^■^^U. See the caption of Fig. |6]for definitions of the 
various plots. 



On the other hand, regarding the average statistic {A'.^{L)), the data from the original NDE as a whole 
are inconsistent: one subset of these data is in good agreement with the GOE model (see Fig. [8]), while the 
complement to this subset is in sharp disagreement (Fig. |9). This inconsistency was not be revealed before, 
as previous analyses ||9][T0l were not undertaken separately for subsets of the NDE. 

To ascertain whether the results shown in Fig. |9]might be due to the inclusion of p-wave resonances, 
we analyzed the data from the corrected NDE after additional exclusion of ^"^^Th and ^■^^U. Results of 
the analysis led to the same conclusion as in the previous case; the validity of the GOE model could be 
rejected with statistical significance of >99.99%. As the excluded p-wave resonances accounted for only 
1.5% of the total number of resonances (of the corrected NDE, excluding ^■^^Th and ■^■^^U), this result was 
not too surprising. 

Results of analyzing the full corrected NDE, are illustrated in Fig. Q As can be seen from this figure, 
exclusion of firm and probable p-wave resonances from the NDE leads to categoric disagreement with 
predictions of the GOE model for L values as small as 70. This disagreement persists even if analysis is 
further restricted to only ^'^^Th and ^asjj resonances of the corrected NDE, with unacceptable disagreement 
again beginning at L « 70. 

It could be argued that the disagreement shown in Fig. |7]is due to missing s-wave resonances and/or 
inclusion of p-wave resonances in the corrected NDE. It is well known that it becomes increasingly difficult 
to correctly determine resonance parity with increasing energy. Therefore, we repeated the analysis of the 
corrected NDE after shortening each sequence by factor of 0.5. In this case we found that at L ft; 72 the 
average quantity {A^^^'^lL)) still deviates appreciably from its expectation value, by « 3.2a. Additional 
Monte Carlo simulations led to the conclusion that within the GOE model the probability (A3(L)) exceeds 
a confidence limit of 3.2cr at one value of L is still very small, 0.04 %. So, even with such a drastic reduction 
the data do not support the GOE model. Acceptable agreement was achieved only at values L < 60. 

5.3 Overall assessment of the NDE 

Given that our analyses of neutron widths as well as energies for the same NDE resonances reveal many 
disagreements with predictions of RMT for the GOE, it is obvious to ask why the full NDE was found 
to be ||9l|T0l in excellent agreement with this theory. As detailed in Ref. HI, the most likely explanation 
has to do with the fact that many of the data in the NDE were selected using this same theory. Almost 
all the neutron data in the NDE were obtained by the group at Columbia University. In their papers, they 
readily admit that they did not have any specific tests for separating s- from p-wave resonances. Hence, 
they routinely used measures derived from RMT for the GOE to perform these separations. Given the 
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numbers of total and assigned s-wave resonances involved, there typically were astronomically large (e.g. 
10^'^ in the case of ^■^^Th ||8]) numbers of possible "s-wave" sets from which to choose. Hence, it does not 
seem surprising that at least one set could be found which agrees with the statistics predicted by RMT for 
the GOE tested in Refs. ||9l[l0|. In other words, the NDE likely demonstrates that it is possible to take an 
incomplete data set and make it agree with GOE statistics by also making it impure. 

With the experiment techniques available at the time, it was just not possible to obtain a complete 
and pure set of s-wave resonances, especially considering that many of the NDE nuclides are very near 
the peaks of the p- and valleys of the s-wave neutron strength functions. Even with the improvements 
in experiments over the intervening 40 years, it still is not possible to obtain such a pure and complete 
set of s-wave neutron resonances as the NDE is purported to be. The main problem remains the reliable 
separation of small (neutron width) s- from large p-wave resonances. In section |8] we briefly describe 
some new experiment techniques that have greatly improved the ability to make firm resonance spin and 
parity assignments. Application of these techniques to nuclides near peaks of the s-wave neutron strength 
function should provide much better data. 

6 Other reported deviations from the PTD 

As noted above, data for the NDE, ^"'^Sm, and i92.i94pj. ^jj g^clude the PTD with high confidence. In 
addition, because the ^"'^Sm and i92a94pj. ^^j.^ demonstrate that can be on either side of the PTD 
value of Vn = 1, the NDE result, which combines analyses of 24 different nuclides, likely underestimates 
the probability with which the PTD can be excluded. There have been a few other reported significant 
deviations from the PTD of which we are aware. 

There have been several reports IITll4l4l that the lowest energy neutron resonances in ^'^^Th deviate 
strongly from the PTD. Our ML analysis ll33l of the data verifies these reports; the data change from 
disagreeing {Vn = 3.8 ± 1.3) to agreeing (z/„ = 0.83 ± 0.68) with the PTD over the first two groups of 
25 resonances. Hence, data for both ^'^^Th and ^''^Sm indicate that the shape of the distribution can 
change in a relatively narrow energy region. As described in Ref. Il34l . our Pt data also suggest that Vn 
may depend on energy. 

The combined data for five odd-A nuclides also were found ifTSl to deviate strongly from the PTD 
despite the fact that very few resonances apparently had been missed. As was the case for ^^^Sm and 
^•^^Th, Vn was found to be significantly greater than 1.0 for these nuclides. 

7 Brief overview of theory 

Our Pt results have inspired a number of theoretical papers [fTGHTSl . some of which have their roots in 
earlier publications (e.g. see Ref. ||201 ). Model calculations ll35l also have shown that deviations from the 
PTD may sometimes occur even when spacing statistics such as A3 do not deviate from GOE predictions. 
However, these latter calculations were for levels near the ground state where collective effects can be 
important and hence leave unexplained how such effects can manifest themselves at the relatively high 
excitations corresponding to typical neutron separation energies. 

In Ref. lfT6l it was proposed that the standard transformation from measured to reduced neutron width, 
r°=r„/\/£^yieV, may be different for nuclides near the peaks of the s-wave neutron strength function. 
However, analysis [36l| of our Pt data showed that the new transformation resulted in increased deviations 
from the PTD. In addition, this proposal could not explain cases such as ^''^Th which are near minima in 
the s-wave neutron strength function. 

Strictly speaking, RMT predictions are valid only for bound states. Possible deviations related to the 
"openness" of the system above the particle separation energy were, to our knowledge, first explored in 
Ref. II20I . Recent models in this vein are described in Refs. 01711181 . The model of Ref. IfTTl does not 
appear to be able to explain our Pt result because the required coupling parameter is ~ 1000 times larger 
than that estimated from our data. The model of Ref. ifTSj results in width distributions broader than the 
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PTD, in agreement with our Pt data. However, it is unclear whether significant deviations from the PTD 
would remain if more realistic matrix elements were used. This model also predicts fluctuations larger than 
the PTD for electromagnetic transitions, in qualitative agreement with our data discussed above. 

As far as we know, no published model can, in general, explain i/„ values greater than 1.0 (although 
the model of Ref. ifTSl can result in f„ > 1 under certain conditions) or large changes in i/„ or within 
relatively narrow energy regions as observed in ^^^Sm and ^■^^Th, and suggested [34] by our Pt data. 

8 New experimental techniques 

New techniques |[3]|22l|25][371 for determining spins and parities of neutron resonances should improve 
the separation of resonances with differing spins and parities and hence help to overcome what has been 
the largest barrier to better tests of the PTD. All these new methods make use of information contained 
in the 7 cascade following neutron capture. It is an old idea ll38l . but improvements in neutron sources, 
detectors, and data acquisition apparatus has led to large improvements over previous results. For example, 
in the case of '^^Mo neutron resonances, applying these methods with a single pair of low-efficiency CeDg 
detectors increased the number of firm resonance assignments from 32 ||39i to 218 Ii25il in this very 
difficult case. 

9 Summary and Conclusions 

We have found significant disagreements with predictions of RMT for the GOE in neutron-resonance data 

main barrier to more and better tests of theory is reliable 
separation of s- from p-wave resonances, but new techniques appear to be well on their way to surmounting 
this barrier For neutron resonances, tests involving the PTD appear to be the most sensitive and reliable, 
largely because experiment limitations can be incorporated into the statistical methods in a straightforward 
manner Tests involving resonance-spacing statistics such as those in Refs. ||9][10| appear to have ignored 
these limitations, and the good agreement between data and theory claimed appears to be due to the fact 
that many of the data used were selected using the theory being tested. Other tests that are purported to be 
less sensitive to missing or spurious resonances have been proposed BOl . but have been shown 01] to lead 
to ambiguous results for sample sizes typical of nuclear data. 
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